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I  INTRODUCTION  AND  SUMMARY 


Evaluation  of  the  performance  of  HF'radio— communication  systems 

requires  that  certain  basic  characteristics  of  both  the  signal  and  the 

noise  with  which  it  is  received  be  known.  Propagation  plays  a  major 

role  in  determining  these  characteristics.  Consequently,  pronounced 

changes  in  HF  propagation  behavior  produced  by  high-altitude  nuclear 

detonations  can  affect  the  performance  of  a  HF  communication  system. 

* 

Complex  computer  codes  have  been  developed  to  assess  the  impact  of  such 
propagation  changes  on  signal  characteristics.  Applying  these  codes  to 
HF  radio  noise  is,  however,  computationally  awkward  because  of  the  wide 
geographic  distribution  and  variability  of  natural  noise  sources.  This 
report  (the  second  of  two^  prepared  for  the  Defense  Nuclear  Agency) 
describes  a  model  that  avoids  this  difficulty  under  conditions  of  suf¬ 
ficiently  severe  and  widespread  propagation  attenuation,  such  as  would 
result  if  several  high-altitude  "blackout"  detonations  were  included  in 
a  nuclear  attack. 

Price  and  Smith^  outlined  a  simplified  HF  noise  model  that  could  be 
applied  to  severely  disturbed  propagation  environments,  and  discussed 
the  approximations  through  which  the  model  was  developed.  A  finished 
implementation  of  the  model  was,  however,  not  presented;  we  complete 
this  work  in  Section  II  of  this  report.  The  data  base  used  to  evaluate 
the  empirical  parameters  through  which  lightning  is  related  in  the  model 
to  equivalent  noise-source  temperature  has  also  been  expanded,  improving 
the  accuracy  of  the  values  obtained  for  these  parameters.  Additionally, 


E.g.,  NUCOM 

^The  first  report  is  Price  and  Smith. ^ 
the  end  of  this  report. 
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the  model  has  been  used  to  calculate  the  effects  on  the  noise  of  a  dis¬ 
turbed  environment,  similar  to  that  found  in  the  vicinity  of  Johnston 
Island  following  the  STAR  FISH  detonation.  Comparison  of  the  model 
results  with  HF  noise  measurements  made  in  Hawaii  shows  satisfactory 
consistency  (the  ionospheric  disturbance  produced  by  STAR  FISH  was  not 
sufficiently  uniform  over  a  wide  area  so  that  the  model  could  be  strictly 
applied  to  these  data). 

A  treatment  of  locally  produced  atmospheric  radio  noise  is  presented 
in  Section  III.  This  noise  contribution,  which  is  unaffected  by  high- 
altitude  detonations,  is  not  properly  dealt  with  in  the  simpler  model 
discussed  in  Section  II.  The  model  departs  from  the  Rayleigh  amplitude- 
distribution  function  (implicitly  assumed  in  the  simpler  model)  at  the 
low-probability/high-amplitude  extreme  of  the  distribution  when  local 
noise  is  included.  These  high  amplitudes  occur  more  frequently  with  the 
improved  model  (although  still  infrequently)  than  is  predicted  by  the 
Rayleigh  distribution. 

The  simple  model  presented  in  Section  II  can  be  readily  incorporated 
into  current  HF  nuclear-effects  propagation-prediction  codes;  such  in¬ 
clusion  is  recommended.  The  local-noise  model  discussed  in  Section  III 
is  not,  however,  suitable  for  incorporation  into  these  codes  in  its 
present  form.  We  recommend  that  this  model  first  be  reduced,  insofar 
as  possible,  to  table  form,  with  the  tables  calculated  using  an  improved 
numerical  version  of  the  present  implementation.  Finally,  extension  of 
the  methodology  developed  to  treat  local  noise  to  provide  a  unified 
treatment  of  both  local  and  remote  noise  is  also  recommended.  Such 
extension,  along  the  lines  outlined  in  Section  III.B,  will  improve  the 
description  provided  by  the  overall  model  of  atmospheric  noise  amplitude 
statistics  in  disturbed  propagation  environments.  Finally,  experiments 
to  validate  the  local-noise  model,  and  its  generalization,  should  be 
undertaken  because  the  experimental  data  upon  which  the  model  is  based 
are,  as  is  noted  in  Section  III  of  this  report,  weak  in  several  particulars. 


II  TESTING,  CHECKING,  AND  EVALUATING  THE  BASIC  NOISE  MODEL 


A.  Introduction 

A  basic  model  for  natural  HF  noise  in  severely  disturbed  propagation 
environments  was  described  in  Price  and  Smith. ^  The  task  reported  on 
here  is  to  review  this  model,  including  the  improvements  that  have  been 
implemented  in  accordance  with  the  recommendations  made  in  Price  and 
Smith. 1  In  this  process,  we  have  accomplished  the  following: 

(a)  Checked  the  physical  parameters  and  confirmed  that 
they  are  correctly  included. 

(b)  Expanded  the  data  base  used  to  derive  the  values  of 
empirical  parameters,  thereby  improving  their  accuracy. 

(c)  Applied  the  model  to  demonstrate  its  performance  in 
its  intended  environment. 

The  physical  parameters  involved  in  the  model  were  indicated  in 
Price  and  Smith, ^  but  simple  averages  and  only  rough  estimates  were  used 
to  establish  the  functional  relationship  between  thunderstorm  activity 
and  equivalent  noise  temperature.  For  ease  in  making  this  relationship 
more  precise,  a  computer  program  representing  the  model  has  been  written 
on  the  SRI  Prime  computer.  Several  changes  from  the  earlier  version  have 
been  included.  The  dependence  of  equivalent  source  temperature  on  eleva¬ 
tion  angle  is  now  included  in  the  lightning  integral,  and  an  analytic 
representation  of  the  normal  D  region  has  been  added. 

The  data  base  for  deriving  the  empirical  relationship  between  thunder¬ 
storm  activity  and  equivalent  noise  temperature  has  been  expanded  to 
include  nine  additional  cases.  A  larger,  more  diverse  data  base  can 
probably  improve  the  model  further,  but  care  must  be  exercised  in  its 
selection  because  only  in  certain  cases  is  the  model  applicable  to  the 
normal  environment.  A  sample  case  in  a  nuclear  environment  was  run  on 
the  program  to  show  how  the  model  will  work. 
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Section  II. B  describes  the  equations  that  are  used  in  the  current 
version  of  the  basic  model.  Section  II. C  describes  the  computation  of 
the  empirical  coefficients  and  evaluates  them.  Section  II. D  describes 
and  evaluates  the  results  of  a  test  case  run  for  a  nuclear  environment 


B.  Model  Equations 

The  basic  model  was  derived  in  Price  and  Smith.'*'  This  section 
explains  in  detail  some  portions  of  the  model  that  were  only  outlined 
in  the  earlier  report.  This  model  is  only  for  severely  disturbed  con¬ 
ditions,  and  is  not  a  general  model.  The  severely  disturbed  environment 
can  be  described  as  one  in  which  D-region  losses  are  high  over  a  large 
region.  This  condition  can  be  produced  in  a  nuclear  environment  if  high- 
altitude  radioactive  debris  is  widespread,  extending  in  all  directions 
for  a  considerable  distance  (thousands  of  kilometers)  from  the  receiver. 
More  precisely,  the  limitations  imposed  are: 

(1)  Only  sources  relatively  near  (within  a  few  thousand 
kilometers  of)  the  receiver  are  significant. 

(2)  Lightning  activity  in  this  source  area  is  (statis¬ 
tically)  the  same  as  that  in  the  immediate  vicinity 
of  the  receiver. 

(3)  Sky-wave  propagation  characteristics  are  uniform 
over  this  region,  with  high  losses;  therefore,  only 
one-hop  paths  are  significant. 

The  general  equation  for  antenna  noise  temperature,  T.(f),  presented 

1  A 
in  Price  and  Smith  (page  13)  can  be  written: 


n/2  2tt 

J  d0  sin  0  ^  d$  p(0,  <p,  f)  £l  +  R2(0,  f)J  Ts(0,f) 


TAf) 


v/2  2it 

/  dQ  sin  0  J 


(1) 


dQ  sin  Q  J  d$  p(0,  4>,  f)  |l  +  R  (0, 
0  0 


[l  +  R2(o,  f)] 
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where 

f  =  operating  frequency  (MHz) 

$  =  azimuthal  angle 
0  =  zenith  angle 

p(0,  p,  f)  =  receiver  antenna  power  gain  pattern 

Tg(0,  f)  =  apparent  sky  temperature 

R(0 ,  f)  =  ground  reflection  loss  (this  parameter 
was  omitted  in  Price  and  Smith. 

The  apparent  sky  temperature,  Tg(0,  f),  is  the  source  term  in  the  noise 
model.  It  is  described  in  terms  of  the  most  important  noise  sources, 
lightning  and  cosmic  processes.  As  a  function  of  zenith  angle,  0, 
these  sources  are  mutually  exclusive.  The  switch  from  one  source  to 
the  other  occurs  as  the  path  of  the  arriving  noise  shifts  from  ionospheric 
reflection  to  penetration.  The  zenith  angle  for  the  switch  is  calculated: 


r  =  radius  of  the  earth 
e 

r  =  distance  from  center  of  earth  to  the  peak  of  the 
p  ionospheric  electron  density  profile 

f  =  critical  frequency  of  the  ionosphere. 

The  galactic  noise  is  received  from  the  region  0  =  0  to  0  ,  and  the 
atmospheric  noise  is  received  in  the  region  0  =  ©c  to  tt/2.  The  apparent 
sky  temperature  from  either  source  is  reduced  by  propagation  losses 
through  the  absorbing  region,  and  the  absorbing  medium  responsible  for 
them  itself  radiates  thermally.  These  phenomena  can  be  represented 
mathematically  as  follows: 


(1)  Sky  Temperature  Due  to  Lightning  Source 

for  0  =  0  (f)  to  7r/ 2 
c 

Ts(0,  f)  =  Tl(0,  f)  [(1  +  R2(0,  f)] 


x 


(f)  sec  i/20 


-2L  (f)  sec  i/20\ 

7^ - ) 


+  T 


-2Lv  (f)  sec  i/20  -2Ly  (f)  sec  i/20 

1  -  I  10  °  +10  X 


where 


L  (f)  =  one-way  vertical  absorption  for  frequency, 
vo  f,  on  ordinary  ray 


L  (f)  =  one-way  vertical  absorption  for  frequency, 
x  f,  on  extraordinary  ray 

i  =  ray  incidence  angle  to  the  ionospheric  layer 
at  the  maximum  absorbing  height  (D  region 
approximately  70  km) 

Tr  =  D-region  radiative  temperature  (283  K) 

(note  change  of  notation  from  Price  and 
Smith  ) 

T  (0,  f)  =  lightning-noise  source  temperature 
-a-, 

=  a„f  L  a,  sec  0 
0  h 

2 

=  lightning  flashes/km  /hour 

=  (0.03  T  +  0.009  Tm2)/720*  *  TA^  *  24/100 
(Cianos  and  Pierce^) 


* 

Hours  per  month. 
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where 


T  =  thunderstorm  days  per  month  (as  given  in 
Reference  3.) 

TA^  =  thunderstorm  activity  per  hour  (as  given 
in  Cianos  and  Pierce  .) 

Oq,  =  empirical  constants  derived  in  next  section. 


(2)  Sky  Temperature  Due  to  Galactic  Sources 
for  0  =  0  to  ©c(f) 


Ts(0,  f)  =  Tc(f) 


-L  (f)  sec  i/20  -L  (f)-  sec  i/20' 

10  °  +  10  X 


+ 


-L  (f)  sec  i/20  -L  (f)  sec  i/20 
v  v 

1  -  1 10  °  +10  X _ 


u) 


where 

Tc(f)  =  100  A2*3  (CCIR4) 

A  =  wavelength  in  meters  =  300/f. 

The  absorption  losses  Ly  (f)  sec  i  and  1^.  (f)  sec  i  are  doubled  for 

o  x 

lightning  to  account  for  the  ascending  and  descending  portions  of  the 
path;  for  cosmic  noise,  these  losses  appear  once,  because  the  ray  path 
only  makes  one  pass  through  the  absorbing  region. 


11 


C.  Calculation  of  the  Coefficients  for  Computing  the 
Sky  Temperature  Due  to  Lightning 

Empirical  coefficients  for  the  sky-temperature  model  were  computed 
by  using  data  taken  in  a  normal  environment,  with  the  requirement  that 
the  data  must  also  satisfy  the  limitations  listed  in  Section  II. B. 

These  conditions  are  satisfied  if  ambient  propagation  losses  are  suf¬ 
ficiently  high.  A  summer  month,  July,  was  chosen  and  appropriate  data 

4 

values  of  FA(f)  from  CCIR  were  found  for  four  separate  cases.  Table  1 
lists  the  data  used.  The  CCIR  data  base  was  collected  between  1957  and 
1961.  An  average  SSN  of  140  was  estimated  for  that  period.  FA(f)  = 

10  log  TA(f)/TQ  where  TQ  =  288  K. 

The  empirical  coefficients  a  =  5.27  x  10^  and  a  =  -7.11  in  the 

0  -a.  1 

functional  relationship  TL(0,  f)  =  ctgf  x  sec  0  were  derived  by  sub¬ 
stituting  the  parameters  from  Table  1  into  Eq.  (1)  and  solving  for 
Oq  and  by  least  squares.  The  relationship  is  linear  on  a  log  scale. 
Figure  1  shows  a  plot  of  the  data  a  =  a^f  ^  as  a  function  of  frequency, 
f,  together  with  the  least  squares  linear  fit  to  the  data  for  the  four 
cases.  The  standard  error  for  this  limited  data  set  is  0.229  on  the 
log  scale. 

A  check  was  made  to  see  if  the  daily  and  hourly  variation  in  thunder¬ 
storm  activity  over  the  four  cases  fitted,  has  contributed  to  the  quality 
of  the  fit.  Elimination  of  from  the  computation  increased  the  standard 
error  of  the  fit  by  14  percent — a  significant  change  given  the  general 
similarity  of  the  cases. 

D.  A  Comparison  of  Model  Results  with  Observed 
Disturbed-Environment  Noise 

The  STAR  FISH  high-altitude  nuclear  detonation  over  Johnston  Island 
on  9  July  1962,  provides  an  interesting  example  with  which  to  compare 
our  model.  Herman ^  reported  that  the  5-MHz  radio  noise  measured  at 
Kekaha,  Hawaii  was  about  12  dB  below  normal  level  at  one  hour  following 
STAR  FISH,  with  a  smaller  decrease  at  higher  frequencies.  Data  also 
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PARAMETERS  USED  IN  DERIVING  THE  EMPIRICAL  RELATIONSHIP  BETWEEN 


Cl an os  and  Pierce' 


are  available  from  Johnston  Island  for  vertical  absorption  and  for  F- 
region  critical  frequency  variations  following  STAR  FISH.  These  data 
have  been  used  as  input  to  our  noise  model  to  obtain  the  expected  noise 
level  at  Johnston  Island  at  one  hour  after  the  shot.  The  results  of  this 
calculation  are  summarized  in  Table  2,  along  with  the  values  of  other 
parameters  that  have  been  input  to  the  calculation. 

If  the  D-region  ionization  produced  by  STAR  FISH  were  sufficiently 
widespread,  the  noise  reduction  calculated  for  Johnston  Island  and  the 
noise  reduction  observed  at  Kekaha  could  be  expected  to  be  similar. 

This,  however,  was  not  the  case;  the  calculated  noise  reduction  at 
Johnston  Island  caused  by  the  detonation  exceeds  that  observed  at  Kekaha 
by  a  considerable  margin.  This  result  is  consistent  with  Herman's^ 
analysis  of  the  noise  data,  in  which  he  finds  thunderstorm  activity  in 
the  Americas,  eastward  of  Kekaha,  to  be  a  major  source  of  the  noise 
received  there.  The  propagation  path  of  this  noise  was  little  affected 
(at  one  hour)  by  the  STAR  FISH  debris,  which,  at  this  time,  remained 
largely  westward  of  Kekaha. 


It 

Unpublished  DNA  reports. 
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Ill  NOISE  OF  LOCAL  ORIGIN 


Noise  of  local  origin,  that  is,  noise  produced  within  line-of-sight 
(or  ground-wave)  range  of  the  receiving  antenna,  is  unaffected  by  the 
changes  in  ionospheric  propagation  loss  with  which  we  have  largely  been 
concerned  in  Section  II.  A  lightning  discharge  at  such  close  distances 
is,  however,  a  relatively  rare  event.  For  example,  a  mean  rate  of  occur¬ 
rence  of  ten  flashes  per  square  kilometer  per  month,  a  value  representa¬ 
tive  of  regions  of  moderately  high  thunderstorm  activity,  implies  only 
about  a  50-percent  probability  of  a  flash  (mean  duration  0.5  s)  being  in 
progress  within  a  250-km  range  (line-of-sight  for  a  source  at  5-km 
altitude  and  a  receiving  antenna  at  ground  level)  at  any  given  time. 

This  sporadicity  of  occurrence  must  be  accounted  for  in  modeling  noise 
of  local  origin. 


A.  Isolated  Flashes 

The  problem  of  describing  locally-produced  noise  can  be  divided  into 
two  parts.  The  first  of  these  is  to  define  the  probability  that  a  local 
lightning  flash  is  in  progress,  and  the  second  is  to  characterize  the 
noise  given  that  a  flash  is,  or  is  not,  in  progress.  Statistically,  if 
p(E)  represents  the  probability  density  distribution  for  the  noise  field 
envelope  amplitude  at  the  receiving  antenna,  then 

p(E)  =  p(E|no  flash)  P(no  flash)  +  p (E | flash)  P(flash)  .  (5) 


In  the  absence  of  a  local  flash,  the  noise  will  be  considered 
Gaussian  within  the  bandwidths  of  interest,  with  a  variance  determined 
by  the  effective  antenna  temperature,  T.  Thus,  we  have 


p(E(no  flash) 


(6) 
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2  4 

where  the  variance,  o^,  is  given  by  CCIR  as 


ol  =  2,82  x  10~34  £2B^-  (V/m)2  (7) 

0 

if  both  the  receiver  bandwidth,  B,  and  frequency,  f,  are  measured  in 
Hertz,  and  the  reference  temperature,  T^,  is  288  K.  This  noise  component 
represents  the  aggregate  of  more  distant  lightning  flashes  or  cosmic 
noise  (or  both)  treated  in  Section  II;  the  effective  antenna  temperature, 

T,  is  defined  as  T^ ( f ) . 

Regarding  the  noise  emitted  by  local  lightning  flashes,  relatively 
few  studies  have  addressed  the  statistical  characteristics  of  the  HF 

g 

radio  emissions  from  individual  flashes.  Oetzel  and  Pierce  characterize 

this  noise  as  a  quasi-continuous  burst  of  radiation  that  continues  over 

the  duration  of  the  flash — typically  some  hundreds  of  milliseconds. 

9 

Horner  offers  a  similar  description  and  summarizes,  in  terms  of  mean 
spectral  amplitude  as  a  function  of  frequency,  various  measurements  of 
source  strength.  These  data  are  reproduced  in  Figure  2. 

Although  significant  fluctuations  in  radiated  power  occur  over  the 
duration  of  a  flash  (most  notably  the  cessation  of  HF  emission  for  some 

Q 

milliseconds  following  a  return  stroke  of  a  flash  to  ground  ) ,  it  seems 
adequate  for  our  purposes  to  treat  the  HF  radiation  as  a  statistically 
uniform  process  throughout  the  flash.  The  noise  also  will  be  considered 
Gaussian,  or  "white,"  over  bandwidths  of  interest,  although  there  is 
some  evidence3^  for  occasional  isolated,  broadband  spikes  of  HF  radiation 
that  may  not  be  well  represented  by  such  a  model. 

The  randomness  of  flash  location,  which  produces  at  the  receiving 
antenna  a  randomness  in  the  mean  power  received  during  a  flash,  must  also  be 

considered.  The  noise  power  incident  on  the  antenna  is  described  by  the 

2  2 
variance,  of  the  noise-amplitude  distribution.  Generally,  is  a 

function  of  the  distance  between  source  and  receiver.  Because  this 

2 

distance  is  a  random  variable,  so  is  o^.  The  net  noise  received  during 
a  flash  is  the  sum  of  the  noise  background  and  that  produced  by  the  flash. 


where 


p(E|flash, 


(10) 


2  2 

The  probability  density  function  for  a  ,  p(a  ),  is  derived  from  that 

,  2 
for  o^: 

p(o2)  -  p(o2  =  o2  -  o2)  .  (11) 


The  times  of  occurrence  of  lightning  flashes  are  taken  to  be 
uncorrelated  and  (statistically)  distributed  uniformly  over  time  in  the 
short  term.  (The  mean  flash  rate  is,  however,  allowed  to  vary  diurnally 
and  seasonally.)  Thus,  the  probability  of  occurrence  of  N  flashes  within 
some  time  interval,  x,  is  given  by  the  Poisson  distribution: 


P(N,t) 


(12) 


where  n  is  the  mean  flash  rate.  If  t  is  set  equal  to  the  mean  duration 
of  a  flash,  then  it  can  be  shown  that  P(N,x)  also  is  the  probability  that 
N  flashes  are  in  progress  simultaneously;  P(0,x)  and  P(1,t)  will  be 
required  for  evaluation  of  Eq.  (5)  from  Eq.  (12). 

P(no  flash)  ■  P(Q,t)  ■  exp(-nx)  ,  (13a) 

P(flash)  =  P(l,x)  *=  nx  exp(-nx)  .  (13b) 

Assume  for  the  moment  that  the  probability  of  more  than  one  local 
flash  in  progress  simultaneously  in  a  region  of  low  thunderstorm  activity 
is  negligible.  Flash  rates  in  regions  of  high  thunderstorm  activity  are, 
however,  sufficiently  great  that  extension  of  the  model  to  treat  over¬ 
lapping  flashes  is  desirable.  Such  extension,  although  analytically 
straightforward,  is  numerically  taxing.  It  will  be  discussed  in  Section 
III.B. 
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To  obtain  the  probability-density  function  for  the  envelope  amplitude 
of  the  incident  noise,  Eq.  (5)  must  be  evaluated.  If  we  insert  into  this 
equation,  the  expressions  from  Eqs.  (6),  (9),  (10),  (11)  and  (13), 
developed  for  its  component  parts,  we  have 


E  I  E 

p(E)  =  exp  (-fit)  — r-expl - ~ 

°0  \  2o0, 


+  nT  exp (-nr) 


Eq,  CIA)  can  be  easily  integrated  with  respect  to  E  to  obtain  the 
probability-distribution  function  for  E: 


P  (E  >  Ex)  «  dE  p  (E) 


exp (-nt )  exp 


00 

exp(-nt)  J  do2  p^o2  =  a2  -  o2j 


2  2 

To  proceed  further,  it  next  is  necessary  to  calculate  p(o1  =  a  - 

2  1 
Oq) .  This  density  function  represents  the  effect  on  incident  noise 

power  at  the  receiving  antenna,  both  of  the  spatial  distribution  of 
flashes  about  the  antenna  and  of  variability  from  flash  to  flash  in  the 
mean  strength  of  the  radiation  from  a  flash.  The  spatial  distribution 
of  radiation  sources  within  a  flash,  although  of  some  interest  (partic¬ 
ularly  the  vertical  extent  because  of  the  impact  of  source  height  on  the 


distance  to  the  line-of-sight  horizon),  will  not  be  addressed.  Present 
11  12 

data  (e.g.  Taylor  and  Rustan  )  are  insufficient,  and  we  cannot  define 
this  distribution  with  any  confidence.  A  "representative"  source  height 
of  5  km  will  be  used  in  the  numerical  examples  that  follow. 

Few  data  have  been  found  from  which  the  variability  in  mean  source 
strength  from  flash  to  flash  could  be  inferred.  However,  in  light  of 
the  large  variation  in  incident  noise  power  associated  with  the  vari¬ 
ability  of  source  location,  this  limitation  is  unlikely  to  prove  serious. 
The  mean  source  strength  averaged  over  all  flashes  is  obtained  from 
Figure  2.  We  assume  that  the  total  radiated  energy  is  spread  uniformly 
over  the  duration  of  a  flash  to  obtain 

a2s  =f-  S2(f)  (16) 

2 

where  a  is  the  variance  of  the  radiated  field  at  unit  distance;  B  is 
s 

the  receiver  bandwidth;  t  the  mean  flash  duration;  and  S(f)  the  noise 
spectral  amplitude  at  unit  distance  for  the  receiver  frequency,  f. 

2  2 

Given  a  ,  a1  is  known  for  a  specific  source  location,  once  the 
s  X 

propagation  loss  has  been  determined.  Generally,  we  have 

°1  “  °s  8^9’  hs’  hP  (17) 

where  the  loss  function,  g,  depends  on  the  geocentric  coordinates, 

(0,  <j>),  of  the  flash  relative  to  the  antenna,  and  on  the  respective 
heights,  h^  and  hf  of  the  source  and  the  antenna  above  the  earth's  sur¬ 
face.  This  relationship  is  assumed  to  be  known. 

The  loss  function  will  also  be  taken  to  include  a  geometrical  factor 

that  specifies  the  loss  associated  with  polarization  and  orientation 

mismatch  between  the  incident  noise  field  and  the  receiving  antenna.  For 

example,  in  the  case  of  a  vertically  polarized  monopole  antenna,  only 

the  vertical  component  of  the  incident  field  is  relevant,  and  g(0,  <(>; 

2 

h  ,  h  )  will  include  the  factor  sin  (9.),  where  0  is  the  zenith  angle 
s  r  A  A 
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of  the  source  (flash)  as  seen  from  the  antenna.  In  terms  of  the  received 
noise  power  at  the  antenna  terminals,  this  mismatch  factor  can  be  recog¬ 
nized  as  simply  the  receiving  antenna  pattern. 


Unfortunately,  the  dependence  of  on  6  and  $  is  sufficiently  com¬ 
plex  for  even  a  simple  propagation  model  so  that  direct  calculation  of 


p(o^)  is  awkward. 


This  difficulty  can,  however,  be  circumvented  by 


using  the  characteristic  function  for  p(o..),  M  2  (iv) ,  to  link  the 
2  1  °1 

distribution  cf  0^  to  that  of  flash  occurrence  in  the  coordinate  variables 
6  and  <p.  This  last  distribution  is  straightforward;  the  flash  density 
for  a  given  range,  6,  is  simply  proportional  to  the  area  at  that  range. 


We  have  then, 


(iv) 


exp  (~ivo 


(18) 


13 

where  the  characteristic  function,  M,  is  defined  as  the  Fourier  trans¬ 
form  of  the  density  function: 


2  2 
p(ox)  exp(ivo1) 


(19) 


Equivalently, 


•fisb 


d$  p(0,  <j>)  exp 


(20) 


2 

where  0^  (0,  <f>)  is  known  by  virtue  of  Eq.  (17). 

2  2  2 

The  probability,  p(o..  =  a  -  aft) ,  needed  in  Eqa.  (14)  and  (15)  fol- 

1  2  U  2  2 

lows  simply  by  substitution  of  0  -  Oq  for  0^  in  Eq.  (18): 
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_1_ 

2tt 


(21) 


00 

/dv 


(iv) 


exp 


M  2  (iv)  is  unaffected  by  this  substitution,  being  calculated  generally 
°1 

in  accordance  with  Eq.  (20). 

The  probability-density  function,  p(0,  4>),  in  Eq.  (2)  is  given  by 


2 

r  sin  6  d@  d 4> 

p (0 ,  <j>)  d0d<f>  »  — -  (22) 

A 

where  A  is  the  area  within  which  the  flash  is  assumed  to  have  occurred 

and  r  is  the  earth’s  radius.  If  attention  is  restricted  to  flashes 
e 

within  a  limited  range,  such  as  those  within  line  of  sight  of  the  receiving 
antenna,  we  obtain 


A  =  2irr^  (1  -  cos  0  )  (23) 

e  m 

where  0  is  the  angular  distance  to  a  flash  on  the  horizon;  for  a  receiving 
m 

antenna  at  ground  height 


cos  0 

m 


r 

e 


r  +  h 


e  s 


(24) 


The  integration  range  for  0  is  correspondingly  restricted  in  the  calcula¬ 
tion  of  M  2  (iv)  from  Eq.  (20);  we  have,  insertings  Eqs.  (22)  and  (23), 


m 


2  77 


M  2  (iv) 
°1 


2tt(1  -  cos  0 

m 


y  J ' d0  sin  0  J " d(J>  exp£ivc^(0,  <f)J 


(25) 
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As  an  example  of  this  procedure,  we  have  calculated  p(o^)  for  the 
vertical  electric  component  of  the  incident  noise  field  at  the  earth's 
surface,  using  a  simple  line-of-sight  propagation  model.  In  this  model, 
the  noise  power  incident  on  the  antenna  from  a  flash  varies  inversely 
as  the  square  of  the  distance,  d,  between  the  flash  and  the  antenna. 
Equation  (17),  then,  becomes,  for  this  model, 


o 


2 

1 


2 

sin  0 , 


(26) 


where  r  is  the  (unit)  distance  to  which  a  is  referenced.  The  source- 
s  s 

receiver  geometry  is  illustrated  in  Figure  3,  from  which  it  can  be  seen 


FIGURE  3  SOURCE-RECEIVER  GEOMETRY 

FOR  LINE-OF-SIGHT  PROPAGATION 
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that  the  source-receiver  distance,  d,  is  related  to  the  angular  distance, 
0,  by 


d2  =  h2  +  2r  (r  +  h  )(1  -  cos  6)  ,  (27) 

S  C  C  b 

where  h  is  the  source  height.  Additionally,  the  antenna  pattern  factor, 
s 

sin  0A,  is  related  to  0  by 


r  +  h 

sin(ir  -  0.)  =  — - - -  sin  0 

a  d 


(28) 


Insertion  of  these  relationships  into  Eq.  (26)  gives 


2  2 

o=o 
1  s 


rs  (re  +  hs)sin  6 


-.2 


h  +  2  r  (r  +  h  )(1  -  cos  0) 

6  6  S 


(29) 


Eq. 


2 

In  this  case,  o^  is  independent  of  4>;  thus,  the  4>  integration  in 
(25)  can  be  simply  performed  to  obtain 


(iv) 


1 

1  -  cos  0 

m 


d0  sin  0 


(30) 


2 

The  variation  of  o1  as  a  function  of  6  given  by  Eq.  (29)  is  shown 

1  2-3 

in  Figure  4  for  a  source  height  h  =  5  km  and  strength  o  *  1,6  *  10 

2  2  s  S 
(V/m)  .  This  value  of  og  corresponds  to  a  spectral  amplitude  S(f)  = 

0.02  Vs/m  tEq.  (16)]  for  a  frequency  f  »  10  MHz  (Figure  2),  a  receiver 

bandwidth  B  =  1  Hz,  and  a  mean  flash  duration  i  =  0.5  s.  The  horizon 

in  this  example  is  at  a  surface  distance  r  0  =  252  km. 

e  m 

2  2  2  2 
The  o^(0)  shown  in  Figure  4  yields  the  p(o^  -  a  -  Oq)  illustrated 

in  Figure  5,  when  used  in  Eq.  (21)  with  a  background  noise  variance 
2  -20  2 

Oq  “  2.8  *  10  (V/m)  ,  corresponding  [Eq.  (4)]  to  a  minimal  background 
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FIGURE  4  VARIATION  WITH  DISTANCE  OF  MEAN  POWER  RADIATED  BY 
A  SINGLE  LIGHTNING  FLASH 


FIGURE  5  PROBABILITY  DENSITY  FUNCTION  FOR  MEAN  RECEIVED  POWER 
IN  A  ONE-Hz  BANDWIDTH  RADIATED  BY  A  SINGLE  FLASH 


2 


temperature  T  =  288  K.  Over  the  range  of  o^  shown  in  Figure  5,  >>  o^, 

2  2  2  2  ^ 
and  consequently  p(o.  =  o  -  o_)  ^  p(a  ).  The  density  function  peaks  at 

2  2  A  u 

a  =0.1  (pV/m)  ,  which  corresponds  to  the  variance  in  the  received  noise 
field  for  a  flash  at  the  most  likely  range,  namely  the  maximum  range  (cf. 
Figure  4).  The  oscillatory  fine  structure  evident  in  Figure  5,  which  leads 
to  some  negative  (and  hence  impossible)  values  for  p  at  the  small  variances, 
is  an  artifact  of  the  Fourier  methods  used  in  the  calculation. 

2 

Use  of  the  p^),  shown  in  Figure  5,  to  calculate  p(E)  and  P(E  >  E^,) 
from  Eqs.  (14)  and  (15)  yields  the  results  shown  in  Figures  6  and  7. 
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FIGURE  6  HIGH-AMPLITUDE  PORTION  OF  PROBABILITY 
DENSITY  FUNCTION  FOR  NOISE-ENVELOPE 
AMPLITUDE  IN  A  ONE-H?  BANDWIDTH  AT 
MODERATE  LtVEL  OF  LOCAL  THUNDERSTORM 
ACTIVITY 


FIGURE  7  HIGH-AMPLITUDE  PORTION  OF  PROBABILITY 

DISTRIBUTION  FUNCTION  "OR  NOISE-ENVELOPE 
AMPLITUDE  IN  A  ONE-Hz  BANDWIDTH  AT 
MODERATE  LEVEL  OF  LOCAL  THUNDERSTORM 
ACTIVITY 
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In  this  calculation,  a  mean  flash  rate  of  1  per  km  per  month,  correspond¬ 
ing  to  n  =  0.08  s-1  within  line-of -sight  range,  was  assumed.  This  flash 

density  is  representative  of  a  region  of  moderate  thunderstorm  activity 

2 

(  5  thunderstorm  days  per  month  ).  Again,  only  that  portion  of  the  total 
range  of  envelope  amplitude  has  been  shown  in  the  figures  for  which  the 
contribution  of  the  noise  produced  by  local  lightning  is  significant. 

This  contribution,  given  by  the  second  terms  of  Eqs.  (14)  and  (15),  is 
approximately  linear  with  n  so  long  as  nx  remains  small  (for  nx  on  the 
order  of  unity  or  greater,  the  effects  of  overlapping  flashes  also  must 
be  considered,  as  described  in  Section  III. 8). 

The  background  contribution  to  the  noise,  given  by  the  first  terms 
of  Eqs.  (14)  and  (15),  lies  predominantely  at  much  smaller  values  of 
E  than  those  produced  by  close  lightning,  giving  a  peak  in  p(E)  at  E  * 
Oq(1.7  *  10"4  pV/m  in  this  example).  The  probability  density  function 
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I 

i 


is  bimodal,  although  the  peak  shown  in  Figure  6  is  much  smaller  in  mag¬ 
nitude  than  the  one  produced  at  smaller  E  by  the  background  noise.  The 
form  of  Eq.  (15)  further  indicates  that  the  integrated  contribution  to 
the  envelope  distribution  from  the  background  noise  exceeds  that  from 
local  flashes  by  a  similar  factor  (i.e.  in  the  ratio  ni:l),  as  also  shown 
by  the  leveling  out  of  P(E  >  E^,)  for  small  (but  >>  aQ)  in  Figure  7. 


B.  Overlapping  Flashes 

Extension  of  the  model  presented  in  Section  III. A  to  include  the 
effects  of  overlapping  flashes  is  not  difficult  in  principle.  We  simply 
generalize  Eq.  (5)  to  a  series  in  which  each  term  corresponds  to  a  given 
number  of  flashes,  N,  simultaneously  in  progress: 


P(E)  P(N)  p(E|N)  .  (31) 

n=0 

The  probability  P(N)  has  already  been  determined  (under  the  assumption 
of  no  correlation  between  flashes);  it  is  given  by  Eq.  (12). 

The  envelope  probability  density  function  given  N  flashes  in  progress, 

p(E|N),  is  defined  just  as  for  N  =  1,  namely  by  Eq.  O) — except  that 

2  2  2  2 
p(E|flash,  a  )  becomes  p(E|N,  a  )  and  p(o  ;  becomes  p(o  |n).  Thus,  we 

have 


■/ 


p(E|N)  =»  /  da  p(e|n, 


o2) 


p(o2|N) 


(32) 


Equation  (10)  also  remains  valid  because  we  are  dealing  (as  regards  the 

field,  not  its  envelope)  with  a  sum  of  Gaussian  random  variables,  which 

2  2 

also  is  Gaussian.  Now,  however,  a  represents  the  sum  of  the  a  for 

2  1 
N  flashes  plus  Oq  for  the  background: 


30 


(33) 


oo 


n=0 


2 

The  variety  of  combinations  of  the  o.  that  can  yield  the  same  value 

2  1 
for  a  is  most  readily  dealt  with  through  the  characteristic  functions 

2  2 
for  the  p(a.),  defined  by  Eq.  (19).  Because  o  is  formed  as  the  sum  of 

1  2 

the  statistically  independent  random  variables  c  ,  its  characteristic 

13  1  2 

function  is  the  product  of  those  for  the  individual  o^: 


N 

m  2  (iv>  =  n  m  2  04) 

C  i=0  ai 

2  2 

With  the  exception  of  Oq,  the  are  all  drawn  from  the  same  population. 
Thus, 


M  2  (iv)  -  M  2  (iv)  M1^  Civ) 
o  oQ  ax 


(35) 


As  previously,  Ma2  (iv)  can  be  readily  calculated  in  terms  of  the  geographic 
distribution  of  flashes,  using  Eq.  (20). 


The  characteristic  function  M  ?  (iv)  for  the  background  component 

°0  2 

of  the  noise  can  be  defined  even  though  Oq  has  a  fixed  value  and  is  not, 

in  the  conventional  sense,  a  random  variable.  If  we  consider  that  a  fixed 
2 

a  results  as  the  limiting  case  when  its  distribution  is  narrowed  to  an 

u  2 
infinitesimal  width,  however,  we  recognize  p(Cg)  as  a  delta  function. 

Thus,  from  Eq.  (19), 


M  2  (iv) 
°0 


6(a^  -  o3)  exp(ivo^) 


exp(ivOg) 


(36a) 


(36b) 
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Fourier  inversion  of  M  „  (iv)  now  produces  p(cr  |n): 


-iv(a2  -  o2)  MN2  (iv)  (37) 

J  °1 

[Equation  (37)  reduces,  as  it  should,  to  Eq.  (21)  for  N  *  1*3 

The  use  of  Eq.  (37)  in  Eq.  (32),  and  this,  in  turn,  with  Eq.  (12) 
in  Eq.  (31)  yields 

p<»  -  ^ 

N=0 


p(o2|N) 


*/■ 


dv  exp 


MNx  (iv) 
Oo 


(38) 


If  the  summation  in  this  expression  is  now  brought  within  the  integrals, 
there  results 


P(E) 


exp(-nt) 

2ir 


00 

/-H 


exp(-E2/2 a2) 


j  dv  expf-iv(c2  -  o2)l  ^ 

—oo  ^  J  M-  n 


I  nrM  7 

L  °i 


(iv) 


n: 


(39) 


N-0 


in  which  the  summation  has  been  reduced  to  the  Taylor  series  representa¬ 
tion  for  the  exponential  function,  whence 


e(E> '  *  f  d°2  7  “*(  £) 

0  . 


00 

t  r  i 

r*  — 

x  /  dv  exp  -iv(a2  -  a2)  exp 

nr  M  (iv)  -  i 

J  L  J 

L  1 

L  1  J 

As  previously,  p(E)  can  readily  be  integrated  with  respect  to  E  to  obtain 
P(E  >  Et): 

00 

p(E  >  ET)  -  -i- y da2  exp^- 
0 

OO 

x  /dv  exp  I- iv  (a  2 


The  above  approach  does  not  depend  strongly  in  its  general  develop¬ 
ment  on  the  particular  assumptions  made  concerning  the  nature  of  the  noise 
radiated  by  an  individual  lightning  flash.  Furutsu  and  Ishida.1^  for 
example,  have  obtained  somewhat  comparable  results  employing  a  model  in 
which  non-uniform  bursts  of  discrete  noise  pulses  are  folded  with  the 
receiver-impulse  response  to  define  the  characteristics  of  the  noise 
received  from  an  individual  flash.  We  have  simplified  this  aspect  of 
the  model  by  assuming  that  the  noise  power  is  distributed  uniformly,  both 
temporally  over  the  duration  of  the  flash,  and  spectrally  over  the  receiver 
pass  band. 

14 

Furutsu  and  Ishida  do  not,  however,  decompose,  as  we  have,  p(E) 
into  an  integration  over  the  joint  distribution  p(E,  o2).  This  step 
88rv*8  isolate  the  influence  on  the  noise  of  propagation  loss,  which 
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we  have  considered  the  dominant  propagation  effect,  and  allows  it  to  be 
dealt  with  in  some  detail.  Lacking  effective  computational  means  for 

lh 

handling  the  variation  of  propagation  loss  with  range,  Furutsu  and  Ishida 
were  constrained  in  application  of  their  results  to  situations  for  which 
this  loss  was  approximately  the  same  for  all  sources. 

Some  preliminary  calculations  have  been  made  employing  Eqs.  (40)  and 
(41),  but  the  results  have  not  been  completely  satisfactory.  The  difficul¬ 
ties  encountered  are  largely  numerical,  and  appear  to  result,  for  the 
most  part,  from  the  embedding  of  the  background  noise  contribution,  whose 
characteristic  function  has  constant  amplitude  out  to  large  values  of  v, 
within  the  Fourier  integrals.  These  difficulties  should  lessen  as  the 
mean  number  of  flashes  in  progress  at  any  time  becomes  large.  This  com¬ 
putational  approach  may  thus  have  application  in  the  joint  treatment  of 
noise  arriving  both  by  ground  wave  and  by  sky  wave  from  a  larger  collec¬ 
tion  region. 

For  the  treatment  of  overlapping  noise  flashes  from  within  line-of- 
sight  range  alone,  the  mean  flash  rate  remains  low,  even  for  highly  active 
thunderstorm  regions;  thus,  direct  summation  of  the  series  in  Eq.  (39) 
is  feasible.  If  this  summation  is  performed  following  Fourier  inversion 
of  the  characteristic  function  for  each  N  on  a  term-by-term  basis,  the 
difficulties  noted  above  are  avoided. 

This  direct-summation  approach  was  employed  for  a  mean  flash  rate 
2 

of  10  per  km  per  month,  corresponding  to  an  area  of  high  thunderstorm 

2 

activity  (20  thunderstorm  days  per  month  ),  to  produce  the  results  shown 
in  Figures  8  through  10.  The  mean  number  of  flashes  in  progress  within 
line-of-sight  range,  is,  in  this  case,  appreciable;  nr  =  0.4.  With  the 
exception  of  the  increase  by  a  factor  of  ten  in  mean  flash  rate,  the 
parameter  values  used  in  this  example  are  the  same  as  those  employed  in 
Section  III. A. 

2 

The  shape  of  the  density  function,  p(o  ),  for  the  net  variance, 
shown  in  Figure  8,  is  little  changed  (even  at  this  flash  density)  from 

that  for  a  single  flash-illustrated  in  Figure  5.  Careful  comparison, 

2 

however,  reveals  a  shift  in  probability  toward  higher  o  values.  This 
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change  can  also  be  found,  on  careful  examination,  in  the  envelope 
probability  density  function,  p(E),  Figure  9,  and  the  distribution  fun¬ 
ction,  Figure  10. 


| 
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FIGURE  8  PROBABILITY  DENSITY  FUNCTION  FOR  MEAN  RECEIVED  POWER 
IN  A  ONE-Hz  BANDWIDTH  AT  HIGH  LEVEL  OF  LOCAL  THUNDER¬ 
STORM  ACTIVITY 
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FIGURE  9  HIGH-AMPLITUDE  PORTION  OF  PROBABILITY  DENSITY 
FUNCTION  FOR  NOISE-ENVELOPE  AMPLITUDE  IN  A 
ONE-Hz  BANDWIDTH  AT  HIGH  LEVEL  OF  LOCAL 
THUNDERSTORM  ACTIVITY 
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FIGURE  10  HIGH-AMPLITUDE  PORTION  OF  PROBABILITY  DISTRI¬ 
BUTION  FUNCTION  FOR  NOISE-ENVELOPE  AMPLITUDE 
IN  A  ONE-Hz  BANDWIDTH  AT  HIGH  LEVEL  OF  LOCAL 
THUNDERSTORM  ACTIVITY 
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ATTN 

J.  Phillips 

ATTN:  M.  Kelly 

ATTN 

T.  Evans 

ATTN 

J.  Newland 

E-Systems, 

Inc 

ATTN 

P.  Komiske 

ATTN 

R.  Berezdivin 

Kaman  Tempo 

Electrospace  Systems,  Inc 

ATTN 

W.  Knapp 

ATTN 

H.  Logston 

ATTN 

T.  Stephens 

ATTN 

DAS I AC 

ESI,  Inc 

ATTN 

W.  McNamara 

ATTN 

J.  Marshall 

ATTN 

R.  Ibaraki 

Linkabit  Corp 

ATTN 

R.  Heckman 

ATTN 

I.  Jacobs 

ATTN 

J.  Lehman 

ATTN 

E.  Tsui 

Litton  Systems,  Inc 

ATTN 

R.  Grasty 

General  Electric  Co 

ATTN:  A.  Harcar 

Lockheed  Missiles  4  Space  Co,  Inc 

ATTN 

W.  Imhof 

General  Electric  Co 

ATTN 

R.  Johnson 

ATTN 

C.  Zierdt 

ATTN 

M.  Walt 

ATTN 

A.  Steinmayer 

Lockheed  Missiles  4  Space  Co,  Inc 

General  Electric  Co 

ATTN 

Dept  60-12 

ATTN 

F.  Relbert 

ATTN 

C.  Old/Oept  68-21 

ATTN 

D.  Churchill /Dept  81-11 

General  Electric  Co 

ATTN 

G.  Millman 

M.I.T.  Lincoln  Lab 

ATTN 

D.  Towle 

General  Research  Corp 

ATTN 

J.  Ise,  Jr 

Magnavox  Govt  4  Indus  Electronics  Co 

ATTN 

J.  Garbarino 

ATTN 

G.  White 
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Martin  Marietta  Corp 

ATTN:  R.  Hefiner 

McDonnell  Douglas  Corp 
ATTN:  H.  Spitzer 
ATTN:  R.  Halprin 
ATTN:  U.  Olson 

Meteor  Communications  Consultants 
ATTN:  R.  Leader 

Mission  Research  Corp 


Riverside  Research  Institute 
ATTN:  V.  Trapani 

Rockwell  International  Corp 
ATTN:  R.  Buckner 

Rockwell  International  Corp 
ATTN:  S.  Quilici 

Santa  Fe  Corp 

ATTN:  D.  Paolucci 


ATTN 

R.  Hendrick 

Science  Applications,  Inc 

ATTN 

Tech  Library 

ATTN 

C.  Smith 

ATTN 

D.  Sappenfield 

ATTN 

D.  Hamlin 

ATTN 

F.  Fajen 

ATTN 

E.  Straker 

ATTN 

S.  Gutsche 

ATTN 

L.  Linson 

ATTN 

R.  Kilb 

ATTN 

R.  Bogusch 

Science  Applications,  Inc 

ATTN 

SZ 

Mitre  Corp 

ATTN 

C.  Callahan 

Science  Applications,  Inc 

ATTN 

B.  Adams 

ATTN 

J.  Cockayne 

ATTN 

A.  Kymmel 

ATTN 

G.  Harding 

SRI  International 

ATTN 

G.  Smith 

Mitre  Corp 

ATTN 

M.  Baron 

ATTN 

M.  Horrocks 

ATTN 

R.  Livingston 

ATTN 

W.  Foster 

ATTN 

C.  Rino 

ATTN 

W.  Hall 

ATTN 

W.  Jaye 

ATTN 

J.  Wheeler 

ATTN 

R.  Leada brand 

ATTN 

W.  Chesnut 

Pacific-Sierra  Research  Corp 

ATTN 

D.  Neil  son 

ATTN 

F.  Thomas 

ATTN 

J.  Petrickes 

ATTN 

E.  Field,  Jr 

ATTN 

R.  Tsunoda 

ATTN 

H.  Brode 

ATTN 

A.  Burns 

4  cy  ATTN 

G.  Price 

Pennsylvania  State  University 

4  cy  ATTN 

V.  Hatfield 

ATTN 

Ionospheric  Research  Lab 

Sylvania  Systems  Group 

Photometries,  Inc 

ATTN 

I .  Kohl  berg 

ATTN 

I.  Kofsky 

ATTN 

J.  Concordia 

ATTN 

R.  Steinhoff 

Physical  Dynamics,  Inc 

ATTN 

E.  Fremouw 

Technology 

International  Corp 

ATTN 

W.  Boquist 

Physical  Research,  Inc 

ATTN 

R.  Deliberis 

Tri-Com,  Inc 

ATTN 

D.  Murray 

RAD  Associates 

ATTN 

W.  Wright 

TRW  Electronics  &  Defense  Sector 

ATTN 

R.  Turco 

ATTN 

R.  Plebuch 

ATTN 

M.  Gantsweg 

ATTN 

D.  Dee 

ATTN 

R.  Lelevier 

ATTN 

W.  Karzas 

Utah  State 

University 

ATTN 

B.  Gabbard 

ATTN 

J.  Dupnik 

ATTN 

H.  Ory 

ATTN 

L.  Jensen 

ATTN 

F.  Gilmore 

ATTN 

K.  Baker 

ATTN 

C.  Grei finger 

ATTN 

P.  .Haas 

Visidyne, 

nc 

ATTN 

J.  Carpenter 

RAD  Associates 

ATTN 

C.  Humphrey 

ATTN 

B.  Yoon 

Rand  Corp 

ATTN 

E.  Bedrozian 

ATTN 

C.  Crain 
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